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Abstract 

We present a brief review of the present day situation with studies of high-temperature superconductivity in iron pnictides and 
'chalcogenides. Recent discovery of superconductivity with > 30 K in A^^Fea A72Se2 (A=K,Cs,Tl,) represents the major new step 
in the development of new concepts in the physics of Fe - based high-temperature superconductors. We compare LDA and ARPES 
data on the band structure and Fermi surfaces of novel superconductors and those of the previously studied FeAs superconductors, 
especially isostructural 122 - superconductors like BaFeaAsa. It appears that electronic structure of new superconductors is rather 
different from that of FeAs 122 - systems. In particular, no nesting properties of electron and hole - like Fermi surfaces is observed, 
casting doubts on most popular theoretical schemes of Cooper pairing for these systems. Doping of novel materials is extremely 
important as a number of topological transitions of Fermi surface near the F point in the Brillouin zone are observed for different 
doping levels. The discovery of Fe vacancies ordering and antiferromagnetic (AFM) ordering at pretty high temperatures (T^ > 
500 K), much exceeding superconducting Tc makes these systems unique antiferromagnetic superconductors with highest T^v ob- 
served up to now. This poses very difficult problems for theoretical understanding of superconductivity. We discuss the role of both 
vacancies and AFM ordering in transformations of band structure and Fermi surfaces, as well as their importance for superconduc- 
tivity. In particular, we show that system remains metallic with unfolded Fermi surfaces quite similar to that in paramagnetic state. 
Superconducting transition temperature T^ of new superconductors is discussed within the general picture of superconductivity in 
multiple band systems. It is demonstrated that both in FeAs - superconductors and in new FeSe - systems the value of T^ correlates 
with the value of the total density of states (DOS) at the Fermi level. 
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1. Introduction 

Discovery of iron based high-temperature superconductors 
HI] attracted a lot of scientific attention leading to a remark- 
able flow of experimental and theoretical works (for review 
see 10, Istl)- The main classes of iron (pnictides and chalco- 
genides) based superconductors known at the moment are: 

1. Doped REllll (RE=La,Ce,Pr,Nd,Sm,Tb,Dy) Fe pnic- 
tides with Tc about 25-55 K, with chemical compositions 
like RE Oi_.vFiFeAs 



IJ. 



2. 



11 1411511^ a 1811210. 

Doped A 122 (A=Ba,Sr), such as Bai_.,KAFe2As2 ll 
Iini and Tc about 38 K. 

1 1 1 systems like Lii_,FeAs with ~ ISKfld,^ 
(Sr,Ca,Eu)FFeAs [.18, 19] with T,. ~ 36 K l^JJ. 
Sr4(Sc,V)206Fe2(P,As)2 with T, ~ 17 K |E~ 



3 
4 
5 

6. FeSe^, FeSei_.,-Te., with T,, up to 14 K Il23l 

7. (K,Cs);cFe2-vSe2 and similar with T^ up to 31K ll24|25[l 



Among these, most recently discovered Fe chalcogenides like 
Kj^Fe2Se2 and CS;cFe2Se2 with rather hi^livalues of supercon- 
ducting transition temperature: Tc=3 IK 124] and 27K [IS] were 
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followed by =3 IK in (Tl,K)FefSe2 126] and form apparently 
a distinct new class. These systems are isostructural to the FeAs 
122 - systems , while Tc values for these Fe chalcogenides 
are a bit smaller than in similar pnictides. According to most 
recent data, in most of the chalcogenides Fe vacancies are in- 
trinsic, so that most general chemical composition is usually 
written now as AAFe2-vSe2. In particular, Fe vacancy ordering 
was discovered in Ko.sFci 6Se2, and most strikingly this system 
seems to be an ordered antiferromagnet (AFM) with pretty large 
Neel temperature (about 578K) iUlH. Further ARPES inves- 
tigations of these systems produced experimental Fermi surface 
maps i29l. fsol fsill significantly different from previously stud- 
ied for Fe pnictides 12, [H. Review of recent findings on Fe 
chalcogenide superconductors can be found in Ref. 11321] . 

In this work we discuss electronic structure, densities of 
states (DOS) and Fermi surfaces of Fe chalcogenides, as com- 
pared to Fe pnictides. We present some estimates on the values 
of Tc, demonstrating that in all Fe - based superconductors there 
is a definite correlation between T^ and the values of DOS at the 
Fermi level. Also we analyze the role of Fe vacancies and AFM 
ordering in KA Fe2-.v/2Se2 in formation of its electronic spectrum 
and Fermi surfaces, which are compared with available ARPES 
data. 
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2. Electronic structure 

There are plenty of papers on LDA band structure of Lai 1 1 
J33, 34, 35], LaOFePjMl, RElll seriesjm BaFezAsz [38, 
[39l|4Q|], LiFeAs ||ill|4i7(Sr,Ca)FFeAs lElEi], Sr42622 Q. 
Electronic structure of Fe(S,Se,Te) materials was discussed in 
Ref. 1 46]. LDA calculations of electronic spectrum of newly 
discovered (K,Cs)j:Fe2Se2 were described recently in ll47[l48ll . 

To illustrate the general picture of the energy spectrum of 
Fe pnictides and chalcogenides in the upper part of Fig. 1 we 
show LDA calculated total, Fe-3d and As-4/7 DOS (left panel) 
matched with the band dispersions (right panel) for typical rep- 
resentatives of both classes. From Fig. 1 one can see that around 
the Fermi level (from -2.5 eV to -1-2.5 eV) there are practically 
only Fe-3d states present, while As-4/7 states are contributing 
at lower energies (from -2.5 eV down to -6.0 eV). 
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Figure 1 : LDA calculated band dispersions and densities of states for typical 
representatives of iron based superconductors: LiFeAs (upper panel) for pnic- 
tides and KFe2Se2 (lower panel) for chalcogenides. The Fermi level Ef is at 
zero energy. 

In a bird eye (large-energy scale) view K vFe2Se2 has similar 
band dispersions as those in pnictides. However, there are some 
quantitative differences, e.g. all Fe-3d and Se-4p states in new 
systems are separated in energy in contrast to Fe-3d and As-4p 



Bal22. Also Se-4p states are of about 0.7 eV lower than As-4p 
states. At the same time, similarly to pnictides the Fermi level 
Ef in Fe chalcogenides is crossed only by Fe-3d states. 

However, at lower energy-scale there is a major difference 
in spectra of both classes. In Fig. 2 we compare LDA cal- 
culated electron spectrum in the immediate vicinity (relevant 
for superconductivity) of the Fermi level for both Bal22 Issll 
system and KfFe2Se2 14811 . To some extent Bal22 bands near 
Ef (upper part of Fig. 2) would match those for KFe2Se2 if 
we shift them down in energy by about 0.2 eV. Main differ- 
ence between old and new systems is seen around F point. For 
KFe2Se2 systems antibonding part Se-4p- band in the Z-F direc- 
tion forms electron-like pocket. In Bal22 corresponding band 
lies about 0.4eV higher and goes much steeper, thus it is quite 
far away from F point. However, if we dope KFe2Se2 systems 
(in a rigid band manner) with holes we obtain bands around F 
point (close to the Fermi level) very similar to those in case of 
Bal22. Namely at 60% hole doping we obtain Fermi surfaces 
with three hole-like cylinders, while stoichiometric KFe2Se2 
has one small electron pocket and larger hole like one near F 
point. Thus, in fact under hole doping we expect several topo- 
logical transitions of the Fermi surfaces |48], which demon- 
strates potentially rich effects of doping this system. 

In general, the electronic structure of new chalcogenide su- 
perconductors close to the Fermi level is significantly different 
from those of FeAs 122 systems. In particular no nesting of 
electron and hole - like Fermi surfaces is observed [48], cast- 
ing doubts on some of the most popular theoretical schemes of 
Cooper pairing developed for iron pnictides. 

3. Anion height and DOS control ? 

It was discovered in Ref. f49^ that superconducting temper- 
ature Tc nonmonotonically depends on anion height Aza with 
respect to Fe layer (see Fig. 3, triangles). Clear maximum is 
seen at about Aza ~1.37A (see also Table 1). Following this 
idea we performed systematic LDA computations of total den- 
sity of states N{Ef) for number of iron based superconductors 
(see Fig. 3, circles and Tab. 1) which have different Aza fs^. 
Here we also add similar results for the new iron chalcogenides 
(K,Cs)Fe2Se2. Nonmonotonous behavior of DOS can be ex- 
plained by hybridization effects. Namely, as a governing struc- 
tural parameter characterizing hybridization strength one can 
chose a-Fe-a angle - an angle between anions (a) and Fe within 
the same tetrahedron. The value of this angle corresponding 
to the strongest hybridization is 109.45°, i.e. for an ideal an- 
ion tetrahedron with Fe in the very center of it. Other crystal 
structure parameters are not changed very much from system to 
system and do not have any transparent dependence of Aza- 

To estimate superconducting critical temperature Tc as a 
function of the Aza one can apply the elementary BCS the- 
ory, where in the expression Tc = I.IAcjds^^^^ corresponding 
N(Ef) enters into the dimensionless pairing interaction constant 
A = gN{Ef)/2 (g is the appropriate dimensional coupling con- 
stant). Taking the Debye frequency a>D=35Q K in rough ac- 
cord with neutron scattering experiments on phonon density of 
states for Lai 11 ifsoll and Bal22 [51] systems, we can find g 
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Table 1: LDA total DOS N(Ef), calculated and experimental T^. for iron based 
superconductors. 
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Figure 2: Top panel - LDA calculated band dispersions in the vicinity of the 
Fermi level for Bal22; Bottom panel - KFe2Se2. The Fermi level is at zero 
energy. 



to fit the experimental value of for Bal22 system since this 
system possesses probably most stable value of (about 38 K) 
with respect to the way of sample preparation and doping. Thus 
we obtain the value of dimensionless coupling constant /l=0.43. 
Then just fixing the value of g as for Bal22 we obtain values 
for all other systems, taking into account only the appropriate 
change of the density of states (Fig. 3, diamonds and Table 1). 
The agreement of these elementary estimates with experimental 
values of T^. is rather striking. 

In fact we do not adhere at the moment to any specific pair- 
ing mechanism, as ojd in BCS expression does not necessarily 
correspond to phonon frequency, it can denote the characteris- 
tic frequency of any kind of Bosons responsible for the pair- 
ing "glue". Our results only show unambiguous correlation 
between the values of superconducting Tc and those of the to- 
tal density of electronic states at the Fermi level for the whole 
class of iron based superconductors (including the new chalco- 
genides), thus supporting the usual BCS-like pairing mecha- 
nism in these systems. The fit using more elaborate expression 
for Tf, like e.g. Allen-Dynes formula, also produces rather sat- 
isfactory results 1.5 2,1 . Relatively high values of effective pair- 
ing couplings, necessary to obtain experimentally observed val- 
ues of Tc, can be understood as due to multiple band electronic 
structure of new superconductors and importance of inter - band 
couplings 115 311 . Special properties of electronic spectrum, such 
as widely popular "nesting" of electron and hole Fermi surfaces 
are not necessary at all. At the same time, inter - band repul- 
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Figure 3: LDA calculated total DOS values N{Ep) (circles, right scale) and 
superconducting transition temperatures Tc (left scale) obtained from simple 
BCS (stars) and experimental values (triangles) versus anion height Aza over 
Fe layer for a number of iron based high temperature superconductors. 



sion of any kind, leading to 5=^ - pairing seems preferable (non- 
phonon) mechanism of pairing, facilitating higher values of Tc 
in both pnictides and chalcogenides. 

4. Vacancies and antiferromagnetism in KFe2Se2 



The recent discovery Il27ll28ll of Fe vacancies ordering and 
antiferromagnetic ordering at pretty high temperature Tai ~ 580 
K, much exceeding superconducting Tc, in the A^Fe2-.v/2Se2 
(A=K,Cs,Tl,...) makes these systems unique antiferromgnetic 
superconductors with highest Tai observed up to now. This 
poses some difficult problems for theoretical understanding of 
superconductivity. Here we discuss the role of both vacan- 
cies and AFM ordering in transformation of band structure and 
Fermi surfaces. 

Ordering of Fe vacancies in the Ko.gFei 6Se2 system lEsIl 
provides VS x VS supercell. Corresponding ordered vacancies 
together with antiferromgnetic order within the Fe-layer 128 1 



3 



e e 




Figure 4: Schematic picture of Kj gFei 5862 Fe-Iayer witli vacancies (rectangle) 
and experimentally observed AFM order (spin up - circles with "+" inside, spin 
down - circles with "-" inside). Translation vectors con'esponding to AFM 
order are ai and 82. Translation vectors for supercell with ordered vacancies 
are b| and b2. 



are presented in Fig. 4. Translation vectors of the supercell be- 
cause of AFM order are ai and ai- Translation vectors corre- 
sponding to supercell with ordered vacancies are bi and bi- 

We performed LSDA calculations for a simplified crystal 
structure of K0.8Fe1.6Se2 given in Ref. ||28||. Our results are 
similar but not identical to other LSDA calculations on this 
system issll ■ For example we obtain metallic AFM state 
with finite density of states at the Fermi level, while energy 
gap forms at lower energies. Because of folding effects auto- 
matically included in LSDA code, it is rather difficult to make 
direct comparison of LSDA band structure for Ko gFci f,Se2 and 
its Fermi surfaces with unfolded LDA bands and Fermi sur- 
faces like those shown above. Also Fermi surface maps from 
ARPES experiments Jl are routinely shown in the un- 

folded Brillouin zone of parent KFe2Se2, like that relevant for 
nonmagnetic LDA 11471 14811 . So these Fermi surfaces are quite 
difficult to compare with extremely folded LSDA ones 154115511. 



To overcome this problem we use simplified (semi)analytic 
model approach first used by us in Ref. ll56tl . To this end we 
fit LDA bands of Ref. [48] with simple parabolic bare spec- 
tra in the vicinity of the Fermi level using the approach of Ref. 
lf56ll and analyze multiple electron scattering (within iron plane) 
by vacancies (CDW) with potential Vi(r) = 2Ai(cosQir -1- 
cosQ2r) and AFM (SDW) potential Viir) = 2A2(cosXir + 
cos X2r). Corresponding vacancy CDW and SDW vectors are 
Qi = 27r(0.4,0.2), Q2 = 27r(-0.2, 0.4), Xi = 27r(0.1, 0.3), 
X2 =2;r(-0.3,0.1). 

Initial (retarded) bare Green function is taken as: 



1 



E - e,(k) + id 



5ij 



(1) 



where /, j are band indices, e,(k) - /-th band model electronic 
spectrum. Because of multiple scattering on vacancies (CDW) 
and AFM (SDW) order, band electron with momentum k could 
be scattered by any of ten possible momenta (See Table 2 of 



Table 2: Table of scattering vectors summation. 
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scattering vectors summation). Namely, electron can preserve 
its momentum, or change it in nine other ways, being scattered 
by Q/; 0/ = -Q/ (Z = 1,2); X,; X, = -X, (Z = 1, 2); Y = (tt, ;r). 
Thus, to find one-band diagonal Green's function G(k, V) = G 
and nine off'-diagonal: G(k ± Q;, k) = Fujy G(k ± X/, k) = <S>nj^ 
and G(k + Y, k) = *F, we end up with system of 10 linear equa- 
tions. Such approach can also be generahzed to the multiple 
band case in a simplified way similar to that used in Ref. 11561] . 
assuming that for both vacancy and SDW scattering both inter- 
and intraband scattering amplitudes are identical. Thus we can 
obtain for the multiple case the system of ten linear equations 
for following Green functions - G'-', F'' <!)'■', , which now 
have two band indices. The equation for diagonal Green func- 



tion IS written as: 



/=1212 



Once we sum up over / we come to 



(2) 



(3) 



/=I2I2 



The rest of other nine equations for G^ - 2, G'^, Fj - 2,- Fj"', 
= Zi ^'/^ = Zi can be obtained using Table 2 of scat- 
tering vectors summation. More details on this approach will be 
presented in Ref. |57], devoted to similar model of hexagonal 
transition metal dichalcogenides, like NbSe2 and TaSe2. 

Finally by solving this system of linear equations we can 
get the diagonal Green's function G'^(k, k) (/, j - band indices) 
and define spectral functions 



A(E,k) = --Im y G"(k,k). 



(4) 



To account for finite experimental spectral resolution we broaden 
all our results by substituting E ^ E + iy, with y - Q.03eV, 
corresponding to typical ARPES resolution. 

The values of model scattering amplitudes Ai and A2 and 
chemical potential were obtained by approximate fitting DOS 
obtained from our model calculations (solid line in Fig. 5) to 
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Figure 5: Comparison of part of total LSDA DOS (dashed line) for 
KvFe2-j/2Se2 with model DOS (solid line) obtained from parabolic bare bands 
with vacancy (CDW) and AFM (SDW) scattering taken into account. Fermi 
level is zero. Letters a,b,c and d denote different doping levels, coiresponding 
to different Fermi surfaces shown in Fig. 7. 
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Figure 6: Model spectral function map for for KvFe2-.v/2Se2 obtained from 
parabolic bare bands with vacancy (CDW) and AFM (SDW ) scattering taken 
into account. Fermi level is zero. 



(similarly broadened) LSDA DOS (dashed line in Fig. 5) in 
a narrow energy interval close to the Fermi level. It can be 
seen that this fit is rather good in this energy interval, where 
parabolic (two-dimensional) fit for energy bands is also satis- 
factory. 

In Fig. 6 model spectral function map (|4]i for KjFe2-.v/2Se2 
is presented. As mentioned above the bare spectra in the vicin- 
ity of the Fermi level were modelled by a number of parabolas 
(cf. Ref. ||56|). The general broadening (finite life time effects) 
comes out from finite y resolution effects. Quasiparticle bands 
with additional shadow like features due to vacancy and AFM 
ordering are clearly seen and can be observed in ARPES ex- 
periments. These bands are to be compared with simple LDA 
results of Ref. [48], shown in the lower part of Fig. 2. 

Model Fermi surface maps for K^Fe2-x/2Se2 with vacancies 
and AFM (SDW) taken into account for different doping levels 
are shown in Fig. 7. Panel (a) of Fig. 7 corresponds to electron 
doped case, panel (b) - undoped and panels (c) and (d) - hole 



Figure 7: Model Fermi surface maps for KvFe2-.i/2Se2 obtained from parabolic 
bare bands with vacancy (CDW) and AFM (SDW) scattering taken into ac- 
count. Panels a,b,c and d correspond to doping levels shown on Fig. 5. 



doped case. Corresponding doping levels are displayed also 
on Fig. 5. These (model) theoretical results are very similar to 
experimental data of Refs. 33, 31 ], with rather large (elec- 
tronic) cylinders around (tt, ;7r)-point and rather small or even 
almost completely smeared out cylinders around F-point were 
observed. Similar picture was found also for parent KFe2Se2 
system in simple LDA [48i]. Presence of vacancy (CDW) and 
AFM (SDW) scattering leads only to formation of some low 
intensive shadow Fermi surfaces which would be rather hard to 
detect in ARPES experiments. In fact, most ARPES investiga- 
tions of these system demonstrate experimental Fermi surface 
maps quite similar to those shown in Fig. 

The most important conclusion is that, rather unexpectedly, 
the role of vacancy scattering and antiferromagnetic ordering 
in the formation of electronic spectrum in the vicinity of the 
Fermi level is rather minor. The system remains metallic with 
well defined Fermi surfaces, allowing the formation of super- 
conducting state on the background of AFM and vacancy order 



5. Conclusion 

In this work we present comparative study of iron based 
pnictide and chalcogenide superconductors. It was shown that 
for both families superconducting transition temperature is ap- 
parently controlled by anion height Az^ with respect to Fe plane, 
which is directly correlated with the values of DOS at the Fermi 
level. Unfolded Fermi surfaces observed in the ARPES experi- 
ments even in strong antiferromagnet K0.8Fe1.6Se2 with ordered 
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Fe vacancies are very similar to those of parent (nonmagnetic) 
KFe2Se2. The system remains metalUc despite AFM and va- 
cancy ordering and superconducting state with rather high 
forms on this unusual background. 
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